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ABSTRACT 

The oceanographic community is moving towards unmanned 
autonomous vehicles to gather data and monitor scientific 
sites. The mission duration of these vehicles is dependent 
primarily on the power consumption of the propulsion system, 
the control system and the sensor packages. 

A customized propulsion thruster is designed. This 
includes a specialized propeller tailored to ABE and a matched 
motor and transmission. A non-linear lumped parameter model 
of the thruster is developed and experimentally verified. The 
model is used to predict thruster performance and compare the 
design thruster with other variants of propeller and 
motor/transmission combinations. 

The results showed that there is a trade-off between 
rapid dynamic response and power conservation. For the 
typical ABE trajectory, the designed thruster provides good 
dynamic response and the lowest power consumption of all the 
modelled thruster units. 
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Chapter 1 



Introduction 



1.1 Motivation 

The Woods Hole Oceanographic Institution (WHOI) is 
currently developing an autonomous underwater vehicle (AUV) 
for scientific survey of the ocean floor. This vehicle has 
been designated ABE, for the Autonomous Benthic Explorer. 

As the oceanographic community explores the ocean floor, 
the reliance on manned submarines to maintain ongoing 
experiments has become restrictive. ALVIN, WHOI ' s manned 
deep submersible, is unable to undertake new research, due to 
an exhaustive schedule maintaining experiments that are 
already in progress. ABE is being developed in order to free 
up assets such as ALVIN and JASON, WHOI ' s unmanned 
submersible, by performing routine data collection and 
surveying at remote ocean bottom scientific sites. 
Additionally, since ABE can operate without a nearby surface 
ship operating as a tender/control ship, ABE will also free up 
valuable research vessel time. In these respects, ABE will 
ccrrplenent the existing capabilities of tethered and manned sufcrrersibles . 
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ABE will be a long endurance vehicle with a typical 
deployment of one year in length. This compares to an on- 
station endurance of eight hours for a manned research 
submersible (ALVIN) and around a month for tethered robots 
(JASON) . During the deployment, ABE will observe a relatively 
small area ( on the order of square kilometers ) of the ocean 
floor at frequent intervals. Figures 1.1 and 1.2 show the 
preliminary ABE configuration and a typical mission profile. 

During a mission, ABE will remain in a semi-dormant state 
for the majority of the time. At a predetermined interval or 
in response to a trigger event, ABE will wake, and preform a 
photographic survey along a preprogrammed flight path. Upon 
completion of the survey, ABE will return to its mooring, 
power down and wait for the next survey time. 
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In order for ABE to have sufficient battery power for a 
one year mission, careful attention must be paid to developing 
a highly efficient propulsion system. This research is 
motivated by the need to develop this propulsion system. 



1.2 Research Objectives 

There are two major components of a AUV propulsion 
system. The first is the mechanical /hydrodynamic system 
commonly referred to as a thruster or propulsor. The thruster 
consists of an electric motor, a transmission and a 
propeller/duct . The second component is the electronics and 
the algorithm used to control the mechanical system. 




Figure 1.2 The Typical ABE mission profile 
will consist of four phases: Descent (1-6), 
Sleep (6), Survey (7-8), and Ascent. [1] 
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The thruster must be optimized for the specific vehicle. 
This involves designing a specialized propeller and duct 
assembly then matching a suitable motor and gearbox to take 
advantage of the motor's high speed efficiency and the 
propeller's low speed efficiency. This is done to maximize 
the conversion of electrical energy to thrust. 

The control system must then be designed to provide the 
best possible dynamic response, while taking advantage of the 
mechanical system's most efficient operating conditions. 

The objective of this research is to develop a propulsion 
system that provides good dynamic response while maintaining 
the high efficiency needed by ABE to perform its mission. 

1.3 Outline of Thesis 

Chapter 2 presents the design of an efficient propeller 
for the ABE vehicle. Chapter 3 examines the selection of a 
motor and gearbox matched to the ABE propeller. In Chapter 4, 
a lumped parameter model of the thruster is developed. 
Chapter 5 contains the experimental verification of the model. 

Chapter 6 compares the steady state performance of the 
designed ABE thruster with other thruster units. Chapter 7 
examines the trade-offs that must be made between power 
efficiency and beneficial dynamics. Chapter 8 summarizes the 
results of the thesis and provides recommendations for further 
research . 
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Chapter 2 



The Propeller and Duct 



The goal of this chapter is to develop a propeller and 
duct combination for ABE that optimize the power conversion of 
the thruster. 

2.1 Introduction to Propeller Design 

The propeller design presented here is based on the 
Massachusetts Institute of Technology's (MIT's) propeller 
development software. This software uses Lifting Line Theory 
and optimum circulation to solve the complex hydrodynamics of 
a marine propeller. [2] 

The process used to design the propeller for ABE is as 
follows : 

1. Calculate the drag of the tentative ABE 
vehicle at the desired operating velocity. 

2. Determine the thrust required. 

3 . Determine the physical constraints of the 
vehicle that effect propeller size. 

4. Estimate the wake field near the propeller. 
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5. Perform a parametric study using MIT's 
Propeller Lifting Line program (MIT-PLL) . 
Minimize Power and torque for a given thrust. 

6. Use MIT's Propeller Blade Design program 
(MIT-PBD10 ) to calculate the blade shape from 
MIT-PLL 's optimum output. 

2.2 ABE Vehicle Constraints 

2.2.1 ABE Physical Considerations 

The ABE vehicle consists of two buoyancy pods supporting 
a single instrument cylinder underneath. The buoyancy pods 
are twenty one inch diameter series 58 bodies, and the 
instrument cylinder is a twelve inch diameter streamline 
cylinder. All three bodies are seven feet long and they are 
connected by a series of struts. There are seven 
thruster/propulsor units on ABE. Three are main propulsors 
mounted at the stern of each cylinder. The remaining four are 
thrusters for attitude and depth control. Two of these last 
four thrusters are mounted vertically, and the other two are 
mounted athwartships . ABE is designed to have a cruising 
speed of one knot, and a minimum sprint capability of two 
knots. The power budget allows 100 watts for normal 
propulsion with a peak of 200 watts available. The propulsors 
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are limited to a diameter of eighteen inches by the space 
available on ABE. 

2.2.2 Vehicle Drag 

The buoyancy pods are made up of twenty-one inch, series 
58, streamlined shapes. Each pod will have a surface area of 
28.1 ft 2 . The drag on this type of body for relatively slow 
speeds is primarily due to friction on the surface. The drag 
can be calculated from the following equation: 

Drag = p U 2 C D S (2.1) 



Where the drag coefficient, C D , is approximately equal to the 
flat plate frictional drag coefficient, C F . For ABE at one 
knot in seawater, the Reynold's Number is =6 X 10 5 and the 
corresponding C F is <0.007. [3] S is the wetted surface area, 
U is the velocity through the water, and p is the density of 
the seawater. Using this equation, the drag of each of the 
buoyancy pods at one knot is 0.56 lbf . The cylindrical 
instrument case will be fitted with streamlined nose and tail 
cones and equation (2.1) holds for this case as well. The 
drag for the twelve inch diameter instrument case at one knot 
is 0.44 lbf. The combined drag of the three main body 
sections of ABE is 1.6 lbf. 
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Attention must be paid to the struts connecting the body 
sections. Careful design of these structures can reduce 
vehicle drag significantly. If it is assumed that the struts 
are well designed and have a streamlined cross-section, then 
equation 2.1 above can be used to calculate the drag. The 
term S is now equal to the entire wetted surface of the strut 
(top and bottom surfaces) . Assuming nine feet of struts with 
one inch maximum thickness and five inches width, the 
calculated drag is =0.15 lbf. 

However, if the struts are poorly designed the drag 
increases by an order of magnitude. Assuming a 1-inch 
circular cross-section, the drag is now calculated by: 

Drag STRUT = ± p U 2 dC d (2.2) 

Where the drag coefficient, C D =1.0, for laminar flow and d is 
the diameter of the strut. For nine feet of circular struts 
the drag is approximately 2.25 lbf. 

For ABE, the struts are assumed to be made up of 
streamlined shapes for two-thirds of the total length and 
bluff sections for the remaining length. The bluff sections 
account for joints, fouling, and imperfections in the strut 
sections. Using this distribution the drag due to the struts 
is 0.8 lbf at 1 knot. 

A large portion of the drag will come from the attitude 
control thrusters. The attitude control thrusters are 



18 



comprised of four thrusters mounted perpendicular to the flow 
streamlines when ABE is traveling in the forward direction. 
These thrusters will be used to make minor corrections to the 
vehicle's depth and heading. These thrusters will be ducted 
to prevent fouling and impact damage. Each unit will have a 
3 inch diameter, streamlined motor case and a 3 inch by 18 
inch duct that will be streamlined for forward motion. Using 
the above equations, each thruster will have a drag of 0.2 
lbf . 

The total vehicle drag for a forward speed of 1 knot is 

3.2 lbf. For the purposes of propeller design, each thruster 
will be designed to give 3 lbf at 1 knot. This allows ABE to 
maintain a full speed capability in the event of the failure 
of one of the three main propulsion thrusters. Additionally, 
ABE will have an excess propulsive force available during 
normal cruising conditions. 

2.3 Optimum Efficiency 

Using the thrust of 3 lbf determined from the drag, the 
vehicle's speed of one knot, an estimated shaft speed of 100 
rpm, and the 1.5 foot propeller diameter, the ideal efficiency 
can be determined from the Kramer diagram, Figure 2.1 [2]. 
The entering arguments for the diagram are X (a form of the 
ship advance coefficient, J s ) and C T , the thrust coefficient. 
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For ABE these coefficients are: 



k = 




Vs 

nn d 



1.6 87 8 — 

sec 

« (100 ipm) (^- — ) (1 . 5 ft) 

6 0 sec 



0.22 



(2.3) 






3 lbf 

4* (1.99 05^2i) (1.6 878— ) 2 ( .75ft) 2 

2 f C * a 



0.60 




[ 2 ] 
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The ideal efficiency for a three blade propeller 
operating under these conditions is rii~83%.[2] This 
represents the maximum efficiency possible, neglecting viscous 
forces . 

2.4 A Parametric Study using MIT-PLL 

Using the information already discussed, we are almost 
ready to begin using MIT-PLL to start developing some chord 
and thickness distributions for the ABE propeller blade. In 
order to enter MIT-PLL, initial chord and thickness 
distributions must be assumed and a inflow velocity field 
needs to be determined. 

The initial chord and thickness distribution were assumed 
to be linear. At the hub, the chord is 2 inches, the 
thickness is Vi inch. At the blade tip, the chord is Vi inch 
with a tenth of an inch thickness. 

For an accurate inflow field to be determined, extensive 
model testing must be performed. Since this type of testing 
is beyond the scope of this paper, a simple inflow field is 
assumed. This field is shown in Figure 2.2. 

With this information, a parametric study of several 
potential propellers was conducted. The parameters used in 
this study are: (1) Number of Blades, (2) RPM, and (3) Blade 
Shape (in a qualitative sense) . The goal is to find the 
maximum efficiency and the minimum power required to generate 
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Figure 2.2 The assumed inflow velocity- 
field 



the required 3 lbf thrust, with an acceptable blade shape. 

The first step is to determine the number of blades for 
the propeller. To do this, the required thrust was varied 
from the 3 lbf value calculated above. This determined two 
things: 1) How dependent efficiency calculations are on 

thrust, and 2) how dependent these calculations are on the 
number of blades. Table 2.1 shows the results of these runs. 
For each run, the optimum RPM was determined, the effective 
wake recalculated, circulation optimized and chord lengths 
optimized. The output was run through the cycle until 
horsepower, torque and efficiency remained constant. Figure 
2.3 . shows the plot of efficiency versus thrust for V s =l knot 
and various numbers of blades. The figure shows that the 
number of blades has little effect on the efficiency over 
ABE's expected operating range. 
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Thrust— » 
Bid: 3 


0.5 

lbf 


1.0 

lbf 


2.0 

lbf 


3.0 

lbf 


4.0 

lbf 


5.0 

lbf 


8.0 

lbf 


10.0 

lbf 


HP 


.004 


.003 


.006 


.019 


.014 


.018 


.033 


.044 


Torque 


.14 


.23 


.39 


.69 


.60 


.82 


1.22 


1.48 


RPM 


45.5 


58.7 


77.9 


92.8 


105 


115 


142 


157 


n 


.875 


.819 


.999 


.702 


.666 


.637 


.573 


.542 


Blades: 4 


HP 


.004 


.009 


.006 


.009 


.019 


. 018 


.033 


.044 


Torque 


.16 


.26 


.39 


.60 


.75 


.91 


1.35 


1.63 


RPM 


41.3 


53.1 


70.1 


83.1 


94.4 


104 


.26 


141 




.875 


.820 


.751 


.704 


.668 


.639 


. 576 


. 544 


Blades: 5 


HP 

i 


.001 


.003 


.006 


.009 


.014 


. 018 


.033 


.044 


Torque 


.16 


.28 


.47 


. 65 


.81 


.97 


1.44 


1.75 


" " ' t 

RPM 


38.9 


48.2 


65.2 


77 .2 


87.7 


96 . 9 


119 


131 


n 


.872 


.818 


.75 


.704 


.669 


.640 


.576 


.545 



Table 2.1 (No Tunnel Used in Calculations) 
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Over the range of thrust from 1 to 6 lbf, the three, four 
and five blade propellers show only a slight variation from 
each other. This difference can be considered statistically 
insignificant. The three, four and five blade propellers were 
used in the next stage of development. All propellers have 
comparable results. 




The evolution of the final blade took many iterations. 
The first 3-blade propeller had a predicted efficiency of 
*70%. (This was calculated without a tunnel around the 
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propeller.) The resulting chord distribution had an extremely 
narrow blade with a very sharp tip. This blade was not 
physically suited for use in the marine environment. Any 
fouling or impact with an obstruction (eg. fish) would have 
destroyed the blade. 

Several chord distributions were used. By taking the 
output of one run, the planer blade shape was plotted (based 
on the chord distribution). From this plot, a new, more 
rugged chord distribution was developed. This chord 
distribution was used for the next MIT-PLL run as the initial 
blade input. For each run of MIT-PLL, the thrust (3 lbf) , the 
diameter (1.5 ft) and the inflow velocity field were held 
constant. The RPM, circulation, and chord lengths were 
optimized. After several iterations, a final 3-blade 
propeller with a 66.1% efficiency and a 4-blade propeller with 
a 66.4% efficiency were chosen. These propellers showed the 
highest efficiency, with durable blade dimensions. 

The three blade propeller is presented since a 
commercially available propeller similar to the designed 
propeller was readily obtained. The commercial propeller is 
presented in section 2.7 . 

2.5 MIT-PBD10, Blade Shaping 

After using MIT-PLL to determine the desired circulation, 
chord and thickness distributions, the actual blade shape and 
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camber that will develop the desired circulation must be 
calculated. The code used for this is MIT's Propeller Blade 
Design program, PBD10. To enter MIT-PBD10, a rake and skew of 
the blade is required in addition to the MIT-PLL output . 
Since rake has negligible effects and serves no purpose for 
the ABE vehicle, no rake is used. The primary purpose of skew 
is to balance (by phase shift) unsteady forces on the 
propeller. Since the forces on ABE are small and the speeds 
of operation are low, the unsteady forces are neglected and a 
small amount of skew (8° at the tip) was added to aid in 
obstruction shedding. 

The recommended default values for MIT-PBD10 were used 
initially. These values determine the nature and extent of 
the wake field. The two dimensional blade cross-section shape 
was chosen as a NACA a=0.8 mean line. Slight modifications 
were made to the extent and contraction of the wake field in 
order for PBD10 to run smoothly in this particular case. The 
PBD10 output fof Kj, Kg and the induced velocities are similar 
to what was described by PLL. Figures 2.4 through 2.7 show 
the resulting blade shape determined by the above process. 
Figures 2.8 and 2.9 show a planer view of the blade, with and 
without skew. 
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The summary of the propeller designed is: 



Number of Blades : 


Z=3 


Hub Radius : 


R h =1 .8" 


Tip Radius : 


O 

II 

CD 


Duct OD: 


a 

n 

CD 


ID: 


d=17 " 


Optimum RPM: 


N=94 .5 


@1 knot & 3 lbf 




Horsepower: 


P= . 008 


Torque : 


t=0 .46 ft-lbf 


Efficiency: 


T|=0 . 661 



Chord Distribution 




Figure 2.4 Chord Distribution 
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Comber (Fo/C) Distribution 




Figure 2.5 Camber Distribution 
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Blade Pitch Distnbuuon 




Figure 2.7 Pitch Distribution 



Planer Blade Shape, No Skew 




Blade Width, Inches 



Figure 2.8 Planer Blade Shape, No Skew. 
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2 . 6 The Duct 

The next step in the ABE propulsor development is to 
design the duct. The primary purpose for the ducts on ABE is 
to prevent fouling of the propellers and the shafts during 
prolonged deployments. The goal is to design a low drag duct 
that has a minimal effect on the vehicle's hydrodynamics and 
propulsion efficiency. The MIT-PLL program traces the neutral 
streamlines at the tip of the propeller. If a duct is placed 
along these streamlines, the duct will have no effect on the 
propeller. If a slight angle of attack is placed on the duct, 
lift can be achieved on the duct section. This lift, when 
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summed around the entire duct, becomes a propulsive force (as 
opposed to drag) . Using MIT-PLL, a duct with a slight 
propulsive force was designed. The neutral nose-tail angle of 
attack is 4.4 degrees. The angle of attack for 5% duct 
propulsion is 3.8 degrees. The cross -sect ional shape for the 
duct is the NACA 0008 Basic thickness form. [4] Figures 2.10 
and 2.11 show this duct. 



31 




Figure 2.10 5% Propulsive Duct 
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2.7 The Experimental Propeller 

The test propeller used in collecting data was a Michigan 
Wheel 18-inch diameter, 16-inch pitch. Sailer ™, 3-Blade 
propeller. This propeller was chosen because it most closely 
matches the chord and pitch distribution of the designed 
propeller. Figures 2.12 through 2.14 describe this propeller. 
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Figure 2.13 Experimental Propeller Pitch 
Distribution 




Figure 2.14 Experimental Propeller Planer Blade Shape 
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Chapter 3 



The Motor and Gearbox 



3.1 ABE Vehicle Constraints 

There are several serious constraints placed on the 
mechanical design of the motor and transmission by ABE'S 
mission environment. The two primary considerations are the 
ambient pressure of the ocean bottom and the limited supply of 
electrical power provided by onboard batteries. 

In order to survive the extreme pressure (10,000 psi), 
the motor assembly must either be pressure tolerant or 
isolated from the pressure. In order to isolate the motor, a 
heavy pressure vessel must be constructed and the motor sealed 
inside. There must be a shaft seal around the output shaft. 
This seal must be leak proof at an extremely large 
differential pressure. Such a seal is expensive and produces 
a large added load on the motor. This load reduces system 
efficiency in a dramatic way. Pressure tolerant motors also 
require shaft seals. However, since the internal pressure of 
the motor casing is maintained at a few pounds above ambient 
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pressure, the differential pressure (DP) is very low. For low 
DPs, the added load is very low. Additionally, seals that 
operate with high differential pressure are prone to slight 
leakage. Any leakage of seawater over an extended period of 
time can result in failure of electric devices. For these 
reasons, shaft seals with high differential pressure should be 
avoided in this application, and a pressure tolerant motor 
should be used. 

The disadvantage of the pressure tolerant motor is an 
increased loss of power due to windage. The windage comes 
from the fluidic drag on the rotor of the motor due to the 
presence of a fill fluid. The fill fluid is a non-conduct ive 
fluid maintained at a few psi above ambient pressure and it 
surrounds and fills the motor. This prevents the highly 
conductive and corrosive seawater from entering the motor. 

Since ABE is a battery powered vehicle with a bus voltage 
of 48 VDC, a DC motor is the obvious choice for the prime 
mover of the thruster. 

The transmission must be selected to match the motor, 
which is most efficient at high speeds, to the propeller, 
which has high efficiency at lower speeds. The gearbox must 
also be of sufficiently high quality and precision to minimize 
the losses due to the gearing. 
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3.2 DC Brushless Motors 



3.2.1 Brushed versus Brushless Motors 

DC motors come in two main types: Brushed and Brushless. 
Brushed motors are the most common DC motor. They have a 
mechanical arrangement of split rings and brushes called a 
commutator. The commutator switches the voltage applied to 
the coils depending on the position of the rotor. This 
switching keeps a positive force on the rotor to keep it 
rotating. A brushless motor relies on electronics to provide 
the proper commutation to the motor based on feedback from an 
external rotor position detector. 

In a high pressure environment, spring loaded brushes 
experience increased wear and a tendency to hydroplane on the 
non-conductive fill fluid. The hydro-planing leads to brush 
chatter and an increased heat load due to the increased 
electrical resistance. Brushed motors have a short life 
expectancy in the high pressure environments. 

Since the commutation on brushless motors is accomplished 
electronically, they do not suffer from any of the above 
problems. They are, however, much more expensive and require 
complicated (and expensive) controllers. Due to the 
importance of longevity and reliability in ABE'S mission 
environment, DC brushless motors will be used. 
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3.2.2 Motor Equations 

In order to evaluate the motor, the efficiency of the 
motor is calculated. The motor chosen for ABE is a Pittman 
elcom © 5100 series DC brushless motor. This motor was chosen 
for size and rated capacity. A schematic of the motor, 
defining the variables used in the equations, is shown in 
Figure 3.1 . 
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Figure 3.1 Motor Schematic 



The motor constants and parameters for this motor are: 



Torque Constant: K T 
Back EMF Constant: K E 
Coil Resistance: R c 



0.173 Nm/Amp 
0.173 V/ (rad/sec) 

4.85 ohms [5] 
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For steady state operation, the following equations describe 
the operation of a DC motor: 

Back EMF : V BACK = K E 

Motor Torque : T m = K T I m 

Motor Current : I m = ( V s - V BACK ) /R c 

Where V s is the supply voltage and co^ is the angular velocity 
of the rotor. Efficiency can be calculated as: 






out 



in 






(3.1) 



From equation 3.1, it can be seen that efficiency 
increases with motor speed. This neglects the effect of 
windage which increase with speed. 



3.3 The Fill Fluid and Windage 

The fluid used for compensation is Halocarbon © 0.8 cSt 
fluid. This is a silicon oil based fluid with a viscosity 20% 
less than water. Of the commercially available, non- 
conductive fluids, this fluid provides pressure compensation 
and has the lowest viscosity. For the purposes of efficiency 
calculation, the fluid flow around the rotor will be 
considered to be laminar flow between two parallel surfaces. 
From testing of other motors with this fluid, a linear 
relationship exists between torque due to drag and angular 
velocity. The proportionality constant in this relationship 
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is K* = 5.175 x 10' 5 Nm/ (rad/sec) . The torque lost (X w ) and the 
power lost (P w ) due to windage are: [6] 



= K w<* u and P w = K w u 2 m 



3 . 4 The Gearbox 

Gearboxes are available in a wide variety of types, gear- 
ratios and efficiency. Gearbox efficiency is highly dependent 
on the manufacturer's tolerances and construction procedures. 
For the Pittman motors, a variety of planetary gearboxes are 
available. These gearboxes will be used to determine the 
desired gear ratio for the ABE thruster. In order to 
calculate the desired gear ratio, descriptive equations of the 
performance of the gearbox must be determined. The equations 
must be written in terms of gear ratio. 

To formulate the equations it is assumed that a planetary 
gearbox is made up of an arbitrarily small stage. A certain 
number of these stages are stacked in order to get the desired 
reduction ratio. Each stage has a specific gear reduction 
(1.1:1) and a specific efficiency (rj t ) . The complete gearbox 
then has a gear reduction ratio of (l.l) n :l and an efficiency 
of Tjj n where n is the number of stages (including fractional 
stages) needed to get the desired reduction. For the Pittman 
gearboxes, two advertised gearboxes are a 4:1 and a 17.33:1 
gear ratio with efficiencies of 80% and 64% respectively. If 
the 4:1 gearbox is used as a baseline, n=14.545 and 1^ = 0.985. 
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If these numbers are used to calculate the efficiency for the 
17.33:1 gearbox, n=29.93 and the efficiency is 63%. Since 
this is a good match, the following equations will be used to 
describe the gearbox: 

N = (1.1)" AND T1 G = (0.985)" (3.3) 

3.5 Matching the Motor with the Propeller 

In order to evaluate a proper match between the propeller 
and the motor, a set of equations describing the motor, 
windage and gearbox must be evaluated for various gear ratios. 
Combining the optimum propeller RPM and torque from chapter 2 
with the component equations developed above, a system can be 
evaluated. Figure 3.2 shows the block diagram of the system. 




Figure 3.2 Thruster System Block Diagram 
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The descriptive equations are: 



T prop = 0.46 ftlbf = 0.62386 Nm 



o D „ D = 94.5 RPM = 9.896 

P zo P aec 



T, = 

1 X G r\ G 



(fl) „ = AT, (j) 

W 1 “ tv G'"prop 



- *1 + ** 






<»m = <*>1 



J = — 

M Tf 

A 



= ^BACK + R-T Z m Y £L4CX " 



N g = (1.1) a Tl c = (.985)” 



JC r - 0.173-^2- K s = 0.173 ———- 

1 tf / rad \ 

V aec ; 



Amp 



iC*. = 5.175X10' 5 ■■ ^ - 

" / rad \ 

V aec ^ 



These equations reduce to: 



= tf s (l.l) n G> prop + JCrJ* 
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■r. = 



p z °p 

(l.l) n iC r (0.985) 



+ (l.l) fl Q 



prop 



with the component efficiencies defined by the following: 
Motor efficiency without windage: 



Power out _ 

Power in V g I m 



Motor efficiency with windage: 



'"Imv 



Power out 
Power in 



T (a) + X 

m * w 

v 3 i m 



Gearbox efficiency: 



tl G = ( 0 . 9 85 ) n 



Total system efficiency: 

n T = g 

These efficiencies were calculated for various gear ratios. 
Figure 3.3 shows a plot of these efficiencies and Figure 3.4 
shows an enlargement of the total efficiency curve. By 
inspection of this figure, the ideal gear ratio is about 10:1. 
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Figure 3.3 System Efficiencies 




Figure 3.4 Total System Efficiency 
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3.6 Sensitivity of the Analysis 

This analysis was carried out using several different 
stage reduction ratios and ef f iciencies . The effects of 
changing these parameters are summarized below. 

1. Changing the stage reduction ratio only from 1.1:1 to 
4:1 had little effect on the gear ratio where the peak 
efficiency occurred. It did effect the height of the peak 
significantly . 

2. Varying the stage efficiency alone (fixed 1.1:1 stage 
reduction ratio) had little effect on the peak location while 
the stage efficiency was above 95%. Below 95% the peak moved 
to lower gear ratios. 

3. Varying both the stage reduction ratio and efficiency 
in a coordinated manner to maintain 80% efficiency at a 4:1 
gear ratio, had no significant effect on both peak location 
and height. 

In general, the model of the gearbox and motor developed 
in this chapter is relatively insensitive to the assumptions 
made about the gearbox. The windage loss dominates the peak 
efficiency curve, provided a gearbox with sufficiently high 
efficiency is used. 



45 



Chapter 4 

The Lumped Parameter Model 
of The Propulsor 



4.1 Introduction 

ABE will incorporate a complex control system designed to 
allow the vehicle to operate independently for periods of up 
to twelve months. In order to accomplish this mission, ABE 
must have extremely efficient thrusters and an efficient 
control algorithm. The control algorithm must be based on a 
simple, yet accurate model of the thruster. In this chapter 
a lumped parameter model of the thruster is developed and 
presented in bond graph notation. A simulation of this model 
was performed using MATLAB, and the results are compared to 
experimental data in Chapter 5. 

4.2 Description of the Thruster Unit 

The thruster under consideration, consists of a DC 
brushless motor, a controller for the motor, a 10:1 reducing 
gearbox, a three blade propeller mounted in a duct. Figure 
4.1 shows a schematic diagram of this assembly. 
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Figure 4.1. Schematic representation of the 
thruster unit. 



4.2.1 The Propeller and Duct Assembly 

The propeller and duct assembly as a unit, has the single 
most significant impact on the overall efficiency of the 
thruster. Therefore, the model of this unit determines the 
accuracy and utility of the overall model. Figure 4.2 shows 
the schematic representation of the propeller and duct. In 
order to describe the hydrodynamic relationships, the 
following simplifying assumptions are made [7]: 

1. The ambient fluid is inviscid, incompressible and of 

constant density. 

2. The gravity effects on the fluid are negligible. 
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3. The flow at the inlet and outlet is parallel to the 
thruster axis. ( irrotational ) . 

4. The only energy storage in the fluid is by kinetic 
energy . 

5. The kinetic energy of the ambient fluid is 
negligible . 

6. The ambient pressure is P 0 and it acts equally at 
the inlet and the outlet. 




In order to model the hydrodynamics, first, consider the 
Kinetic Co-energy, E K * . This is the typical "Physics Book 
Kinetic Energy" { Vi mv 2 } . 

= < 4 -d 

Where pV is the mass of the fluid enclosed by the duct. 
However, this volume must be corrected to take account of the 
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added mass effect of the surrounding fluid. This correction 
involves a concept known as Added Mass or Virtual Mass . Added 
mass is a phenomenon that occurs when a body (or a fluid 
element) moves through a fluid. An additional inertia effect 
is added to account for the effort required to "push" the 
fluid out of the way of the moving body. For the case of a 
body moving slowly through a static fluid, the added mass is 
equal to the volume of the body times the density of the fluid 
(ie. the mass of the displaced fluid) [8] and [9]. For this 
model (where the density is constant) this is done by setting 
V to twice the actual enclosed volume. Q is the volumetric 
flow rate. A is the cross-sectional area of the duct. Now, 
define the pressure momentum, T, as: 

r & -j-(Eg) = - P (4.2) 

d£? ' K A 2 

***Note that there is a linear relationship between T and 
Q. This is analogous to the standard translational definition 
of momentum, p=mv, and leads to the numerical equality of the 
Kinetic Energy, E K and the Kinetic Co-energy: 



e- = Jr do = -|pv(-§) 2 



(4.3) 









a 2 t 2 

2p V 



A 2 p VQ 

2 p V A 2 
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In Summary: 



p V&_ _ A 2 T 2 
2 A 2 2 p V 



(4.5) 



r = ■£_££ And 




(4.6) 



Now consider a power balance for the propeller and duct. 
Power in: 

The power flow into the duct comes from three major 
components: 1) The driving motor, 2) Any Kinetic Energy 

flowing into the thruster inlet, and 3) The velocity/ 
opposing-force product at the inlet ( work done by the fluid 
entering) . 

1. The power input from the motor/gearbox is the product 
of the torque, t, and the angular velocity, co. 

2. Since it is assumed that the ambient fluid is at rest, 
there is no kinetic energy flow into the thruster. 

3. The velocity/ force power is the product of ambient 
pressure P 0 and the cross-sectional area and the fluid 
velocity at the inlet. 
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(4.7) 



P in = TO) + 







Power out : 

The power flows out of the duct by two processes. The 
first is the force/velocity work at the discharge. The second 
is the kinetic energy of the discharged fluid. This latter 
flow is called the convected kinetic energy. 

1. Since the cross-sectional area is constant and the 
average velocity of the fluid is the same at the inlet and 
outlet, this term is the same as the corresponding term at the 
inlet, P 0 Q. 

2. The convected kinetic energy is the kinetic energy per 
unit volume times the volumetric flow rate: 



Note: the absolute value preserves the sign of the convected 
kinetic energy to allow for flow reversal. 



K = Convected Kinetic Energy = — |C| = — - 1 ^1 (4.8) 

V |V| 2 pV 2 




(4.9) 
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